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ABSTRACT 

We investigate the star formation properties of ~ 800 sources detected in one of the 
deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~ 0.1 — 4) 
and about four orders of magnitude in star formation rate (SFR). It includes both 
star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into 
radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared 
luminosity, as traced by Herschel , with the SFR computed from their radio emission. 
We find that the radio power is a good SFR tracer not only for pure SFGs but also in 
the host galaxies of RQ AGNs, with no significant deviation with redshift or specific 
SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs 
population and the occurrence of AGNs in sources with different level of star formation. 
Finally we discuss the possibility of using deep radio survey as a tool to study the 
cosmic star formation history. 
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1 INTRODUCTION 


The faint radio sky is a complex mixture of star form¬ 
ing galaxies (SFG) and active galactic nuc l ei (AGN) (e.g., 


Mauch & Sadler 20Q3; ISmolcic et al.l 120081 : IPadovani et al.1 


20091 . 201 ll) . Indeed, radio emission can either be due to 


the relativistic jets powered by the AGN or to synchrotron 
emission from electrons accelerated by supernova explo¬ 
sions. In the latter case, radio emission c an then b e used 
as a star formation r ate (SFR) tracer (e.g., iYun et ahll200ll : 
IPannella et al.l 120091) which opens the possibility of using 
radio survey to study the cosmic star formation history 
(CSFH). Moreover, the radio band has the advantage, com¬ 
pared to UV or optical frequencies, to offer a dust unbiased 
estimate of the SFR being almost unaffected by dust extinc¬ 
tion. 


However, the majority of radio surveys available up to 
now were only sensitive to the most powerful radio loud 
(RL) AGNs, whose radio emission is mainly due to the 
relativistic jets rather than to star formation. In the few 
cases in which the sensitivity was high enough to detect 
SFGs, only the most extreme starburst galaxies, with SFRs 
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> 1, 000 M 0 yr _1 , could be detected above redshift about 1, 
at the mjy level. It is only going below 0.1 mjy, that we be¬ 
come sensitive, over a wider redshift range, to the bulk of the 
SFG population, which moreover represents the dominant 
contribution to th e overall radio population at these radio 
flux densities (e.g., Padovani et all201ll : iBonzini et al .iRoTR : 
IPadovani et al.ll2015l) . Such a sensitivity has been reached so 
far only on small patches of the sky as with the VLA sur¬ 
vey of the Extended Chandra Deep Field South (E-CDFS) 
considered in this work. 


A further problem that has limited so far the use of 
radio surveys for studying the CSFH is the challenge of sep¬ 
arating the two radio emission mechanisms, jets and star 
formation, in faint radio samples. In particular, the origin of 
radio emission in the so-called radio-quiet (RQ) AGNs has 
been a matter of debate for quite some time. It has been 
proposed that t hese sources are a scaled down version of RL 
objects (e.g.. iMiller et all 1 1 99.il : iGiroletti fe Panessal 120091) 
or that the radio emission is mostly due to the star forma¬ 
tion in the host galaxvle.g. JSopp fe Alexanderlll99ll ). In the 
former case, the SFR derived from the radio power would be 
overestimated compared to the one obtained by other SFR 
tracers like, for example, the far-infrared (FIR) luminosity. 


If instead radio emission in RQ AGNs is due to SF, 
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the radio power could be used as an SFR tracer even when 
the AGN emission strongly contaminates the optical-to-MIR. 
host galaxy light. This would allow us to include powerful 
AGN hosts in SFR studies therefore superseding the lim¬ 
itation of many current works that are restricted to pure 
star forming systems. However, it is first necessary to find 
an effective method to separate RQ and RL AGNs, as in 
the latter source s the radio lumi nosity primarily traces the 
jets emission. In iBonzini et al.l (120131) we have proposed a 
method to separate these two AGN populations. We note 
that the RL AGNs contribution to the cosmic star formation 
history is expected to be low for two main reasons; firstly, 
RL AGNs represent only a small fraction (10%) of the radio 
sources at flux densities of ~ 30 pjy and, secondly, their 
host galaxies have usually low star formation rates as they 
are mainly passive systems. 

Moreover, estimating the SF activity in AGNs host 
galaxies is particularly important to investigate the possible 
impact of the AGN in shaping its host galaxy properties. In¬ 
deed, both from the theoretical and observational side, there 
has been a lot of debate on the role the AGN can have in 
regula ting the SF activity, e i ther triggering or suppressing it 


Ji PI 

(e.g.. Mullane^etaL 2012 1 _ Rosario_etjd. 2012; IPaee et al.l 


I2012L lHarrison et alJ 20121 : Zubovas et al.l 20131 ). 

In this work we exploit the star-formation properties 
of a large sample of radio sources selected in a deep VLA 
survey of the E-CDFS down to a 5<r flux density limit of 
about 32 /ijy. The paper is organized as follows. In section 
[5] we described our radio sample, how it splits in the dif¬ 
ferent source populations, and the ancillary data used for 
the analysis. After describing the method adopted to esti¬ 
mate the FIR luminosity of the radio sources (Sec. 0, we 
present the radio-FIR correlation for radio selected SFGs 
(Sec. 2J. In Section [u] we compute the SFR from the radio 
and FIR luminosity and we compare them in Sec. [2] The po¬ 
sition of our radio selected sources in the SFR-stellar mass 
plane is described in Sec. [7] and in Sec. [H] their SF activity 
is investigated. In Sec. 151 and HOI we discuss and summarize 
our results, respectively. Finally, Appendix[A]is dedicated to 
the description of the empirical model we compare with our 
observations and in Appendix m we describe the catalogue 
with the physical properties of the radio sources used in this 
work that we make publicly available. 

In this paper we assume a cosmology with Ho = 70 
km s" 1 Mpc -1 , II m = 0.27, and Da = 0.73 and a Chabrier 
initial mass function (IMF: IChabriedf2003h . 


2 SAMPLE DESCRIPTION 

2.1 Radio data 

The Extended Chandra Deep Field South (E-CDFS) has 
been observed with the Very Large A rray ( VLA) at 1. 4 
GHz between June and September 2007 llMiller et al.ll2008lh 
The survey reaches a best rms sensitivity of 6 pjy and the 
average 5 a flux density limit is 37 pjy with near-uniform 
coverage. A catalogue including sources down to peak flux 
density of five times the local rms noise has been extracted. 
A description of the s urvey strateg y and the data reduc¬ 
tion details are given in I Miller et al.l (l2013l h The radio cata¬ 
logue includes 883 sources. Using a likelihood ratio technique 


dCiliegi et al.ll2003l) . the optical/in frared counterparts of the 
radio sources have been identified (IBonzini et al.l 20121'). Ex¬ 
cludin g the outermost region as defined in IBonzini et all 
(12012th the wealth of multi-wavelength data available for the 
E-CDFS allows rob u st estimates of ph o tometric redshifts 
llSantini et al] I200S : iTavlor et aD 120091 : ICardamone et all 
l20ld : Rafferty et al l201llh The area with photometric red- 
shifts coverage includes 779 radio sources that will consti¬ 
tute our main sample for the analysis presented in this pa¬ 
per. Combining the photometric redshifts with the output of 
several spectroscopic campaign in the E-CDFS, we were able 
to assign a redshift to a to tal of 675 radi o sou rces, 37% of 
which are spectroscopic (see I Bonzini et afll2012l. for details). 
The average redshift is ( z ) ~ 1.1. 


2.2 Radio source populations 


As mentioned in Section [T] the sub-mjy radio population 
is a mixture of SFG and AGN. The latter further divides 
into RQ and RL AGN. To separate these different popu¬ 
lations we use a multi-wavelength approach combining ra¬ 
dio, mid-infrared, and X-ray data. A detailed descriptio n 
of our classification scheme is given in iBonzini et al.) d2013h ; 
here we briefly summarize its main characteristics. We select 
RL AGN using the 52406 s parameter, which is the logarithm 
of the ratio between the observed 24 /im flu x density and 
the observed 1.4 GHz flux density fe.g.. ISargent et aDl2010l . 
and references therein). In the 52406 s-redshift plane we have 
defined an “SFG locus” below which we find radio sources 
with a radio excess compared to the typical 52406 s ratio of 
SF systems. This excess is the signature of an AGN con¬ 
tribution to the radio luminosity; t hese sources have bee n 
referred to as radio-domi nant AGNs (lAppleton et alJl2 004L 
radio-excess AGNs (e.g. iDrake et al.l 20031: De^^oroet_aL| 
2013h or simply radi o-loud AGNs fe.g. iDonlev et al.ll2005l : 


Padovani et al. 120091 4 in the literature. Here, we adopt the 


nomenclature RL AGNs. The SFG locus has been obtained 
computing the 524063 as a f u nction of red shift for the star 
forming galaxy M82 ([Polletta ii etaUj200’7]) and assuming a 
0.7 dex dispersion. In lBonzini et all d2013h we have discussed 
the implications of choosing the M82 template for source 
classification (see also [B]). Within and ab ove this locus , both 
SFGs and RQ AGNs can be found (e.g. IDonlev et al.ll2005l : 
IBonzini et alJl2013h . A source is classified as RQ AGN if 
there is any evidence of AGN activity in the other bands con¬ 
sidered: we classify a radio source as RQ AGN if a) it has a 
hard band [2-10 keV] X-ray luminosity greater than 10 42 erg 
s _1 ; b) it lies in the “AGN wedge” o f the IRAC color-color 
diagram, as defined by IDonlev et al.1 d2012l l. Otherwise, the 
object is classified as a SFG. According to this scheme, our 
sample of 779 radio sources includes 167 RL AGNs, 188 RQ 
AGNs, and 424 SFGs. 


2.3 Herschel data 

The E-CDFS has been observed by Herschel as part of the 
PEP (PACS Evolutionary Probe) programme llLutz et ali 
l201llL The Herschel observations are deeper in the central 
part of the field, the GOODS field, and shallower in the 
outskirts. The whole field has been observed at 100 and 160 
pm and for the central part 70pm data are also available. 
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For GOODS proper we use the com b ined reduction of PEP 
and GOO DS- Herschel llElbaz et alJl201 lh PACS data as de¬ 
scribed in iMagnelli et alJ 1 20131 ). The 5a flux level of the 
100 (160) pm maps are 0.85 (2.1) mjy and 6.25 (13.05) 
mjy in the GOODS and E-CDFS, respectively. The 70 
pm 5a flux density limit is 1.35 mjy. Herschel photome¬ 
try was performed through point spread function (PSF) fit¬ 
ting, adopting Spitzer MIPS 24 /im detected sources as po¬ 
sitional priors (IBerta et al.ll2013h . Blind catalogues are also 
extrac ted by means of PSF -fitting using the StarFinder IDL 
code (Di olaiti et al.l 2000a. b) and include all sources with 
S/N > 3a. We cross-correlate this Herschel/ PACS cata¬ 
logues with our radio sample using a 1.5" searching radius. 
We find a total of 490 matches, 33% in the GOODS field 
and 67% in the outer region of the E-CDFS. 

In particular, in the PACS-radio sub-sample we have 42 
RL AGN, 130 RQ AGN, and 311 SFG. We note that in the 
PACS sample we recover 69% of the R.Q AGN, 78% of the 
SFG but only 25% of the RL AGN from the original radio 
selected sample considered in this paper. The small fraction 
of RL AGN is not surprising since they are preferentially 
hosted in passive, dust poor, galaxies. We therefore expect 
only a small fraction of them to be detected in the FIR. 

For the radio sources without a PACS counterpart we 
compute 5 a upper limits from the local rms noise, therefore 
taking into account the non uniform coverage of the Herschel 
observations in this field. 


We chose this library mostly for its wide range of SED 
shapes. This is important if we want to fit the photometry of 
our radio selected sample since, as already mentioned, it is 
highly non-homogeneous: it includes both normal star form¬ 
ing and starburst galaxies, objects with the clear presence 
of a powerful AGN and sources whose SED is dominated 
at all wav elengt hs by the emission from stars. However, as 
the IBerta et al.l J2013h library is based on PACS detected 
sources, it lacks galaxies with very old stellar population 
and consequently very low emission in the FIR. Therefore, 
we decided to add three templ ates of el liptic al gala xies from 
the SWIRE template library (IPolletta et a! .1 [20071 ). 

We use a standard \ 2 minimization technique to find 
the best fit template. Examples of the optical-to-FIR SED 
fitting for the three class of sources are shown in Fig. [T] 
For sources without Spitzer and/or Herschel detections we 
imposed the best fit template to not exceed the 5 a upper 
limits in these bands. We finally compute the total FIR lu¬ 
minosity (Lpir) integrating the best fit template between 8 
and 1000/rm. Errors are computed repeating the fitting pro¬ 
cedure 100 times for each source, randomly modifying the 
photometry within the uncertainties and then taking the 
standard deviation of the Lfir distribution obtained. 

We note that for sources without PACS detection the 
FIR luminosity is derived from the best fit based only on 
the optical-MIR photometry, taking into account the upper 
limits in the FIR. These luminosity estimates are therefore 
less robust than the one constrained by the PACS measure¬ 
ments. 


3 FIR LUMINOSITY OF RADIO SOURCES 


The Herschel photometry is crucial to have a good esti¬ 
mate of the FIR emission since it allows to trace the cold 
dust emission. At a typical redshift of ~ 1, the PACS mea¬ 
surements probe the rising part and/or the peak of the rest- 
frame FIR. dust emission bump. We apply a fitting technique 
to the full UV-to-FIR spectral energy distribution (SED) 
to better constrain the FIR luminosity of our radio sample, 
taking advantage of the exquisite multi-wavelength coverage 
available in the E-CDFS. 

The UV- NIR photometry i s obt ained combining the 
BVR selected ICarda mone et ahl (l2010h catalo gue, the K se¬ 
lected iTavlor et ah ( 20091 1 catalogue, and the lPamen et all 
(1201 lh IRAC selected SIMPLE catalogue. In the MIR, 
we used the MIPS 24 and 70/rm observations from the 
Far- Infrared Deep Extrag a lactic Legacy Survey (FIDEL) 
llDickins on fe FIDEL Te airJl2007il . A likelihood ratio tech¬ 
nique ( Ciliegi et al.l 120031 ) followed by a careful visual in¬ 
spection has been used to reconstruct the optical-to -MIR. 
photometry of the radio sources (see iBonzini et al.l I2013l . 
for details). At longer wavelengths, given the lower surface 
density of sources, we simply cross-correlate the Herschel 
and radio counterparts catalogues, as described in the pre¬ 
vious section. In the GOODS field, the PACS 70/im detec¬ 
tion rather than the MIPS one is adopted when available. 
In summary, for each radio source we have on average 16 
photometric points. 

To fit the UV-to-FIR photometry we use the 
IBerta et al.l d2013h template library that is based on the ob¬ 
served SED of PACS detected sources. It is a collection of 
32 templates, mostly SFGs, of which twelve also include an 
AGN component. 


4 RADIO-FIR CORRELATION 

The radio and FIR luminosity in SFGs follow a tight 
empirical rel ation, the so-called "radio-FIR correlation" 
(RFC) (e .g., iDickev fe Salpeterl 1 19841 : Ide Jong et al.l Il985l : 
iBdi 1 200311 . 

The physical reason for its tightness and linear¬ 
ity over several decades in luminosity has been in vesti ¬ 
gated by several theoreti cal models (e.g., Lacki__gUal. 1 120101 : 
lLacki fe Thompsonl l20ld : ISchleicher fe Beckll2013ll . but re¬ 
mains still unclear. These works also predict that the RFC 
should break at high z where inverse Compton losses start 
to dominate over the non-thermal synchrotron emission. 

From the observational point of view, there are some 
hints of possible deviations at high z or in extreme ob¬ 
jects like sub-millimetre galaxies, nevertheless these obser- 
vations are still consistent with no sign i ficant evolution (e.g. , 


lvison et al 2Qlol: 


MoricetaL | l20ld : ISargent et alJ 120101 : 


Magnelli et aljhoiS : iLut j|2014l . and reference therein). 


Deep radio surveys would be needed to investigate the 
behaviour of the RFC at high redshift for normal SFGs. 
Indeed, even to detect star forming galaxies with SFR of 
hundreds of solar masses per year at z>2, a pjy radio sensi¬ 
tivity is required. Such a sensitivity has been reached so far 
only for small patches of the sky like in the VLA survey of 
the E-CDFS considered in this work. 

In this section we therefore investigate the radio-FIR 
correlation for our radio selected SFGs up to z~ 3. In Fig- 
ure[2]we show the radio power versus the FIR luminosity for 
SFGs with a clear detection in at least one PACS filter since 
their FIR luminosity estimate are more robust. With the ex- 
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Figure 1. Examples of the optical-to-FIR SED fitting for the three class of sources. The combined photometry is plotted as red symbols 
as a function of wavelength (observed frame) and the solid line show the best fit template. Source ID and redshift are given in the bottom 
right corner. 


ception of two clear outlier^], the two luminosities are tightly 
corre lated and lie along the empirical RFC (e.g., Kcnnicutt 
11998T) plotted as dot-dashed line and parametrized by: 

log(Pl.4GiTz) = log(Z/FIR.) + 11.47 (1) 

where Piaghz is the radio power in W Hz -1 at 1.4 GHz 
and Lfir is the FIR. luminosity in unit of solar luminosity 
expressed in erg s -1 . The average dispersion is 0.2 dex. We 
note that such a dispersion is smaller compared to the dis¬ 
persion in q 24 obs values for our SFGs that is 0.33 dex. This 
implies that the correlation between the total FIR luminos¬ 
ity and the radio power is stronger that the one obtained 
using only MIR data. The color scale in Fig. [2] represents 
the redshift of the sources. Thanks to the /rjy sensitivity of 
the VLA observation, we can detect SFGs in a wide redshift 
range, from z^O.l to z~ 3 and we find that the RFC holds 
over the whole redshift range with almost constant disper¬ 
sion. Moreover, our data probe four orders of magnitudes in 
luminosity and therefore include both normal SFG as well 
as the most active systems with SFR of a thousand of solar 
masses per year. The SF in these latter systems is thought to 
be triggered by episodic violent events such as major merg- 
ers rather th an by secular processes as in normal SFGs (e.g., 
IDaddi et alj|2010bj ). Our study suggests that the correlation 
between radio and FIR emission is the same in both kinds 
of systems, at least down to the luminosities probed by our 
observations. However, we still cannot exclude a breakdown 
at lower star formation su rface densities as predicted by e.g. 
iSchleicher fe Beckl J2013I ) for which even deeper radio data 
would be needed. 

This results justify the use of both the FIR and radio 
luminosity as independent but equivalent SFR tracers for 
our radio selected SFGs. It will be a useful element in our 


1 RID 521 (dark point) is a resolved edge-on spiral whose 24/rm 
emission has been probably underestimated due to aperture pho¬ 
tometry. As a consequence the fit is poor and the FIR luminosity 
is underestimated. RID 577 (light point): there is a small offset 
between the radio and optical-MIR emission hence it is possible 
that the radio emission is associated with a background object. 
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Figure 2. Radio power at 1.4 GHz versus FIR. luminosity for 
SFGs detected by Herschel. The color scale rep resen t the r edshift 
of the sources and the dashed line shows the iKennicu 3 (1199811 
relation given in Eq. |T] 


investigation on the origin of the radio emission in RQ AGN 
(see Sec. ED>. 


5 STAR FORMATION RATE ESTIMATE 

The FIR luminosity is a tracer of recent SF activity as the 
UV radiation from young stars is absorbed by the dust and 
re-emitted at longer wavelengths, in the FIR. The Lfir pro¬ 
vides a better estimate of the SFR compared to the one de¬ 
rived from the UV or line emission since these measurements 
are highly affected by dust extinction. Moreover, comput¬ 
ing the SFR. from the FIR is particularly useful for bright 
AGN hosts. Indeed, the AGN emission can overshine the 
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emission of the whole galaxy at optical wavelength, while 
it represents a tiny fraction (<10%) of the tota l lum inosity 


in the FIR in most of the cases (e.g.. 

latziminaoglou et al. 

2010j: Rosario et al. 2012: Berta et al 

2013f). We adopt the 

Kermicutt (1998|) relation to compute 

:he SFR from the FIR 

luminosity: 


SFR Fm [M e yr- 1 } = 4.5 x lQ-^Lmnlergs- 1 ]. (2) 


As this equation assumes a Salpeter IMF, we rescaled it for 
a Chabrier IMF multiplying by a factor of 0.6. 

Also the radio continuum emission can be used as a 
tracer of recent star formation in SFGs since it is nearly 
all due to synchrotron emis sion from relat ivistic electrons 
associated to SN remnants (I Condon! 11992tl . The empirical 
conversion between the ra dio power at 1 .4 G Hz and the SFR 
of the galaxy according to lYun et alj |200ll ) is: 

SFR^Mq yr- 1 } = (5.9 ± 1.8) x 10 -22 Pi. 4G H*[W.ff2 -1 ].(3) 

We multiply the derived SFR by a factor of 0.6 to convert 
from a Salpeter to a Chabrier IMF. 


6 RADIO EMISSION IN RQ AGN 

In Figure [3] we compare the SFR computed from the FIR 
luminosity with the SFR derived from the radio luminosity 
for all the 675 radio sources with redshift. Different colors 
represent different classes of objects: green for SFGs, blue 
for RQ AGNs, and red for RL AGNs. Sources detected in 
at least one Herschel band are plotted with full symbols. 
We first concentrate on these sources since their SFRpm is 
more reliable. 

As described in Section [2 we find that the SFGs follow 
the RFC over the whole range in luminosities accessible to 
our survey. As a consequence, the two SFR estimates are in 
agreement over four decades in SFR with a typical disper¬ 
sion of 0.2 dex. More interestingly, we find good agreement 
between the SFRs derived from the two different tracers also 
for RQ AGN with only a slightly larger scatter of 0.23 dex. 
The results of a linear regression fit (least square bisector) 
are: 

log(SFR Fm ) = 0.97 ± 0.02 x log (SFR,) + 0.02 ± 0.04 (4) 
and 

\og(SFR Fm ) = 0.96 ± 0.05 x log (SFR. T ) - 0.09 ± 0.10 (5) 

for SFGs and RQ AGNs with PACS detections, respectively. 
The two relations are fully consistent within the uncertain¬ 
ties. This implies that the radio power is a good tracer for 
the SFR not only for SFGs, but also for RQ AGNs, with 
almost the same degree of accuracy for the two types of 
sources. We emphasize that this result is not a direct conse¬ 
quence of our selection method; indeed, even if both types 
of sources lie within the “SFGs locus” in the 924 obs-redshift 
plane they have different MIR characteristics. SFG, for ex¬ 
ample, have an average dispersion of q 24 obs values of only 
0.33 dex while for RQ AGNs the dispersion is two times 
larger (0.65 dex). This is mainly due to the large AGN 
contribution at MIR wavelengths in many RQ AGNs. As 
a consequence, an estimate of the SFR from the observed 
24ym emission would lead to a much more uncertain re¬ 
sult for RQ AGNs as for SFGs. Instead, based on Fig. [3] 


we can claim that the radio luminosity is as good tracer as 
the FIR of the star-formation in RQ AGN host galaxies. 
The correlation shown in Fig. [2] also strongly suggests that 
the main contribution to the radio emission in RQ AGN 
is therefore due to the SF in the host galaxy ra ther than 
black hole activity (see also IPadovani et al .11201 il l. As men¬ 
tioned above, that could not be claimed considering only 
MIR data as these wavelengths are not reliable for deriving 
the SFR in many powerful AGNs; it would have been im¬ 
possible to resolve if the larger scatter in q 24 obs values for 
RQ AGNs compared to SFG was due to the presence of hot 
dust heated by the AGN or to the presence of small radio 
jets. Hence, it was crucial to include Herschel data in our 
analysis that allowed us to obtain a reliable estimate of the 
SFR in our sources, assuming that the AGN contribution to 
the total F IR luminosity is in the large majority of the cases 
negligib le (lHat ziminaoglou et all l20ld : iRosario et al.l 120121 : 
IBerta et al l 1 201 dT see also below). 

The behaviour of the RL AGNs further supports our 
claim that the origin of radio emission in the two types of 
AGN is different: they scatter out from the one-to-one re¬ 
lation in Fig. [3] since the SFR computed from the radio lu¬ 
minosity is overestimated due to the jet contribution to the 
radio power. Note that this is true not only for very pow¬ 
erful RL objects but also for low-power radio AGNs. The 
comparison of the two SFR tracers allows us also to isolate 
sources that have been misclassified; indeed there are some 
RQ AGNs (~ 5%) with PACS detectio ns that are b elow the 
relation, especially at high redshift. In lBonzini et ali l|2013l h 
we estimated a contamination to the RQ AGNs population 
from RL AGNs with a strong contribution from the AGN to 
the 24 ym flux density that boosted their q 24 obs value into 
the SFGs locus of about 5%, in agreement with our findings. 
Therefore, we believe that most of the outliers are indeed RL 
AGNfl These outliers are also the main responsible for the 
slightly larger scatter (0.23 dex) of RQ AGNs as compared 
to the SFGs population (0.2 dex). 

6.1 AGN contribution in the FIR 

As already mentioned in Sec.|3j twelve of the SED templates 
in the library used to fit the photometry of our radio sam¬ 
ple include an AGN component. The AGN contribution to 
the FI R luminosity in these templates has been estimated 
in iBerta et alj d2013l ) . We have therefore the opportunity to 
subtract the AGN contribution to the SFR.p ir and analyse 
its impact on our conclusions. About 50% of the RQ AGNs 
and 35% of the SFGs are best fitted by a template that 
has an AGN component. However, it is important to note 
that the fit solutions are highly degenerate; in particular, 
templates with small AGN contribution and pure SFG ones 
might provide equivalent good fits to the source photome¬ 
try. Hence, the results reported in the following are valid 
only in a statistical sense. Nevertheless, the fact that only 
half of our RQ AGNs are best fitted by an AGN template 
shows that a classification based only on the optical-to-FIR 
properties can miss a significant fraction of AGNs. For these 

2 We will keep the original classification in the rest of the paper. 
Nevertheless, we checked if our conclusions change reclassifying 
these sources, finding no significant differences. 


© 2015 RAS, MNRAS 000, [T] 7? 


































6 Bonzini et al. 



Figure 3. SFR derived from the FIR luminosity versus the SFR 
from the radio power at 1.4 GHz. SFGs are plotted as green cir¬ 
cles, RQ AGNs as blue squares and RL AGNs as red triangles. 
Full symbols represent sources detected in at least one PACS 
filter, while sources shown as empty symbols are Herschel non¬ 
detections. Large symbols with lighter colors are the results of the 
stacking analysis (see Section |6.2| for details). 


Figure 4. Ratio of the SFR computed from the fit to the whole 
photometry and the one obtained using only the UV-to-24 pm 
data as a function of redshift. Different colors correspond to dif¬ 
ferent classes of sources as in Fig. [3] Large yellow diamonds rep¬ 
resent the median and the standard deviation in different redshift 
bins. 


latter sources, the AGN contribution to the optical-to-FIR 
SED is typically minor and, in most of the cases (80%), they 
are identified as AGNs due to their high X-ray luminosity. 
On the other hand, the fact that we find some SFGs that are 
best fitted by an AGN template does not mean that they 
are all real AGNs. As already highlighted, the SED fitting 
technique often yields degenerate solutions, especially in the 
case of discrete libraries. However, it is possible that some of 
these SFGs are indeed AGNs since our classification scheme 
was aimed to select a clean sample of RQ AGNs and a con¬ 
tamination from low- luminosity AGNs in SFG population is 
expected (see I Bonzini et al .1120131 . for details). 

For all the sources fitted with an AGN template we sub¬ 
tracted the AGN contribution to the FIR luminosity and 
re-computed the corresponding SFR. The new best linear 
regressions have a slope of 0.95 for both SFGs (±0.02) and 
RQ AGNs (±0.05) and the offset from the 1:1 correlation are 
0.04±0.05 and —0.10±0.11, respectively. They are therefore 
consistent with the previous estimates within the uncertain¬ 
ties. The only difference is a slightly larger scatter around 
the best fit relation of 0.21 dex and 0.25 dex for SFGs and 
RQ AGNs, respectively. 

We conclude that the AGN contribution to the FIR 
emission as computed from the best fit template has negli¬ 
gible impact on our results. We will therefore consider the 
SFRfui non-corrected for AGN contribution in the rest of 
this paper. 


6.2 Stacking of Herschel undetected sources 

As already discussed above, for sources without a PACS de¬ 
tection the estimate of the FIR luminosity is less robust. 
Indeed, the best fit template has been obtained considering 
the photometry only up to 24pm, with upper limits at longer 
wavelengths. To estimate how this can impact the SFR es¬ 
timate, we have re-computed, with the same method, the 
SFR for sources with PACS detection but considering only 
their photometry up to 24 pm. In Fig. [5] we plot the ratio 
of the two SFR measurements, considering or not the PACS 
data, as a function of redshift. Large yellow symbols rep¬ 
resent the median and the standard deviation in different 
redshift bins. We observe a large number of outliers (~25%) 
and a slight tendency to underestimate the FIR emission 
especially above z^l (~ 0.16 dex) when the Herschel data 
are not included in the fitting procedure. 

This test explains the larger scatter of the empty sym¬ 
bols in Fig|!I]and justifies their exclusion from the linear re¬ 
gression analysis. However, it could be that the good agree¬ 
ment we find between the two SFR tracers holds only for 
the FIR-brightest objects while it breaks for the Herschel 
undetected sources. 

Therefore, we pe rformed a mean stacking analysis 
llBethermin et al.ll2010l) of the PACS-undetected sources to 
estimate their average SFRyib,. We divide the undetected 
sources in classes and for each class we split the sample into 
two bins of radio power (21.5 < P r < 23.5 and 23.5 < P T < 
24.5). We then stacked both the 100 and 160pm Herschel 
maps at the position o f the PA CS-u ndetected radio sources 
using the error maps llLutz et al.ll201ll l as weights. Flux den¬ 
sities are extracted from the stacked images through PSF- 
fitting; aperture corrections and correlated noise corrections 
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Table 1. Results of the stacking analysis. 


(1) 

Class 

(2) 

# sources 

(3) 

mean z 

(4) 

mean P r (i. 4G Hz) 
[log (W Hz-1)] 

(5) 

/l00/i,m 

mJy 

(6) 

/l60/xm 

mJy 

(7) 

S/Nioo^m 

(8) 

S / N 160/xm 

SFG 

21 

0.86 

23.12 

7.78 

12.90 

2.9 

2.7 

SFG 

41 

1.88 

23.98 

2.39 

8.08 

3.6 

2.9 

RQ AGN 

12 

0.68 

22.90 

1.49 

6.48 

1.7 

4.2 

RQ AGN 

24 

2.03 

24.04 

4.88 

11.10 

2.9 

2.4 

RL AGN 

49 

0.57 

22.90 

1.66 

3.04 

2.6 

1.9 

RL AGN 

37 

1.25 

23.98 

1.09 

2.85 

1.3 

1.4 


are applied l|Berta et al .11201.'ill. Uncertainties are e stimated 
using a bootstrap approach llBethermin et al.l20ld l and cor¬ 
rected also for high-pass filtering effects. Tabic [T] summarize 
the stacking analysis results. We obtain > 2cr detections for 
SFG and RQ AGN in both bands in each bin, with the only 
exception of the 100/rm stacked image for RQ AGN in the 
lower radio power bin. We note that this bin has also the 
smallest number of sources (12) and therefore the S/N ra¬ 
tio is lower. RL AGN are detected only at 100/im in the low 
radio power bin while they are undetected (< 2er) elsewhere. 

To compute the average FIR luminosity for each sub¬ 
sample of sources, we build their average optical-to-FIR SED 
and then we proceed as described in Sec. [3] Stacked flux 
densities below the 2cr threshold are treated as upper limits 
in the fitting procedure. 

The corresponding average SFRfib are plotted in Fig[jJ] 
as large full symbols, with colors slightly lighter than the 
ones of the their corresponding source class. The aver¬ 
age SFRsfib. of PACS-undetected SFGs and RQ AGNs lie 
within the scatter of the corresponding linear relations. We 
can therefore conclude that the two SFR tracers are equiv¬ 
alent for our populations of RQ AGNs and SFGs regardless 
of their IR. brightness. 

The average SFRfib of the RL AGNs not detected by 
Herschel (the majority of the population) has been derived 
taking into account the upper limits obtained with the stack¬ 
ing technique in the PACS bands. Those average values are 
significantly smaller than the average SFR. r obtained us¬ 
ing the mean radio power in each bin. This confirms that 
their radio luminosity has a significant contribution from 
the jets and therefore cannot be used to estimate their star- 
formation rate. We also note that, in both low and high 
power sources, RL AGNs have SFRfib. lower compared to 
the other two classes of sources meaning that the host galax¬ 
ies of RL AGNs are on average more passive. 


7 SFR VERSUS STELLAR MASS 

The S F activity of a galaxy can occur in two d ifferent modes 
(e.g., iDaddi et alJ l2010al : iGenzel et al.l l20icil l: a starburst 
one, probably triggered by major mergers or in dense SF 
regions; and a quiescent one, associated with secular pro¬ 
cesses, which is observed in the majority of the SFGs. In this 
second mode, the SFR is correlated with the stellar mass of 
the galaxy, forming the so called “main sequence” (MS) of 
SFGs. The MS has nearly the same slope both at low and 
high redshift but the normalization increases o f ab o ut a fac- 
tor of 20 from the local Universe to z « 2 ('e.g. jNoeske et all 


2007 

; Daddi et al. 2007|; Elbaz et al. 2007|; Pannella et al. 

05 

O 

O 

CM 

: Peng et al. 201(i; Rodighiero et al 20 1 ll: Wuvts et al. 

2011 

)• 


In this framework, it is then interesting to locate out¬ 
sources on the SFR-stellar mass plane. In Figure [5] we plot 
the SFRfib as a function of the galaxy stellar mass for our 
objects divided in four redshift bins. Having already veri¬ 
fied that the radio and the FIR luminosities are equivalently 
good tracers of the SFR for our SFGs and RQ AGNs, we de¬ 
cided to use in the following the SFR. derived from the FIR 
luminosity because it allows us to carry on the analysis also 
for RL AGNs for which the SFR r would instead be mean¬ 
ingless. Moreover, since the SFRfib is most widely used in 
the literature, it makes easier the comparisons with previous 
studies. Stellar masses are computed using a two component 
(AGN+ga laxy) optical-t o-MIR. SED fitting technique as de¬ 
scribed in lBonzini et al.l (120131) . assuming a Chabrier IMF. 

The horizontal dotted line in each pane0 of Fig. [5]marks 
the minimum SFR ( SFRmi n ) accessible to our survey, i.e. 
the SFR that correspond to the flux density limit of the VLA 
observations at the minimum redshift in the specified bin, 
based on equation [3] and corrected for a Chabrier IMF. Note 
that SFGs and RQ AGNs mostly lie above these lines, being 
radio emission in both classes associated to star formation 
related process. RL AGN instead are detected even when 
they have SFR well below the SFRmin since their radio 
emission is boosted by the presence of the jets. 

The dashed lines indicates the position of the MS at 
the average redshift of the sources in each bin. We adopt 
the following law for the redshift evolution of the MS: 

log (SFR(M, z )) = -7.77 + 0.79 log(M) + 2.8 log(l + z) (6) 


where M is the stellar mass expressed in unit of solar masses. 
T he slope and the redshift evolution are based on the results 
of iRodighiero et al.l (1201 if) and the normalization provides 
the best agreement with our data (see Appendix El for de¬ 
tails) and it is c onsistent wit h the normalizations adopted 
in the litera t ure (Noeske et al. (12007 1 for th e local universe; 
lElbaz et al.l d2007h and Daddi et al] d2007h for z ~1 and 
z ~2, respectively). The dot-dashed lines above and below 
the MS correspond to ±0.6 dex. 

Due to our radio flux density limit, we are able to detect 
MS galaxies with low stellar mass (M < 1O 1O ' 5 M0) only in 
the first redshift bin while at higher redshift we probe the 
bulk of the main sequence only at the high mass end. 

We note that RQ AGNs occupy the same locus in the 


3 In the first redshift bin the minimum SFR is 0.07 Mq yr 1 
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Figure 5. SFR versus stellar mass for different redshift bins. The symbols are the same as in Figure[3] The horizontal dotted lines marks 
the minimum SFR accessible to our survey in each redshift bin (see Text for details). The dashed lines indicates the position of the MS 
at the average redshift of the sources in each panel. The loci of starburst and passive galaxies are delimited by the dot-dashed lines (±0.6 
dex) above and below the MS respectively. 


M s tar — SFR plane as SFGs suggesting that the majority of 
the host galaxies of radio selected RQ AGNs are not signif¬ 
icantly different from the inactive galaxies population (see 
Sec. 18.31 for a more detailed discussion). The fact that the 
relative fraction of SFGs over the RQ AGN decreases as a 
function of redshift is due to a combination of selection and 
evoluti onary effects as furthe r discussed in the next session 
and in lPadovani et al] (120151 ). 


8 SPECIFIC STAR FORMATION RATE 

Since the SFR is correlated with the stellar mass, a use¬ 
ful quantity to describe the SF regime of a galaxy is its 
specific SFR (sSFR), i.e. the SFR divided by the stellar 
mass. Sources with sSFR particularly higher with respect 
to the MS are undergoing extremely intense star forma¬ 
tion activity, possibly triggered by major mergers, while 
passive galaxies are characterized by very low sSFRs. In 


this work, we will define as "starburst" those sources whose 
distance with respect to the MS (i.e., A \og(sSFR)MS = 
(og[sSFR(galaxy)/sSFRMs{M 3t ar, z)]) is larger than 0.6 
dex, and we will call “passive galax ies” those with 
A log (sSFR)ms < 0.6 dex (similarly to llRodighiero et all 

Era)- 

In the following, we will discuss the star-formation prop¬ 
erties of the three classes of sources introduced in Sec. E3 
We will treat RL AGN separately as their properties are sig¬ 
nificantly different from the other two classes and the flux 
density limit of the VLA survey has a less direct impact on 
the minimum star-formation rate detectable in these objects 
(see Sec. Q • 

8.1 SF in RL AGNs 

RL AGNs are generally though to be hosted in massive red 
and dead galaxies. Indeed, the majority (56%) of our RL 
AGN are hosted in passive galaxies but a significant frac- 
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Redshift Redshift Redshift Redshift 


Figure 6. Fraction of starbust (magenta circles), main-sequence galaxies (green diamonds), and passive galaxies (red squares) versus 
redshift. The panels refer from left to right to RL AGNs, R.Q AGNs, SFGs and mock catalogue of star forming sources. In this latter 
panel empty (light colors) symbols refer to the mass selected sample, while full (dark colors) symbols to the mock sample "observed" 
with the same radio flux density limit of our VLA survey. 


RL AGN RQAGN SFG Mock catalog 



Figure 7. Fraction of starbust (magenta), MS galaxies (green), and passive galaxies (red) versus stellar mass. The panels refer from left 
to right to RL AGNs, RQ AGNs, SFGs and mock catalogue of star forming sources. In this latter panel empty (light colors) symbols 
refer to the mass selected sample, while full (dark colors) symbols to the mock sample "observed" with the same radio flux density limit 
of our VLA survey. 


tion of them are in MS star-forming galaxies and a small 
fraction (5%) is even hosted in starburst systems. Examples 
of powerful radio loud AGNs hosted in starb urst galaxies 
have been found at high redshift (z>3) (e.g., Ilvison et al.l 
[20121) . We have investigated whether the fraction of actively 
star forming host galaxies increases with redshift. As shown 
in the left panel of Fig. [6] we find no significant trend up to 
z~3. However, we see a strong trend with the stellar mass 
of the host (see left panel of Fig0. Passive galaxies are the 
vast majority of RL AGN hosts at masses > 10 10 ' 5 Mq while 
at lower masses they are mostly MS galaxies. This trend is 
similar to what is observed in the overall galaxy population; 
the fraction of active galaxies in c reases as the stell a r mas s 
decreases (e.g., IPeng et al.l l20ld : lllbert et al.l 120101 . l2013h . 
Therefore, it seems that there is no strong connection be¬ 
tween the presence of radio jets and the SF activity of the 
host. On the other hand, we observe a large spread in stel¬ 
lar masses at low radio power, while the most powerful R.L 
AGNs are only found in the most massive objects (see Fig. 
El suggesting a link between the maximal energy that could 
be released by the jet and the mass of the black hole, assum¬ 
ing the Mbh — M bulge relation. That could also explain why 
in shallower radio surveys the fraction of RL AGN hosted in 
actively star forming systems is generally lower than what 
we observe. As our sources are selected in the radio, it is 
very unlikely that we are missing bright radio sources with 
small stellar mass. 

Finally, it is also unlikely that we are simply misclassi- 



21 22 23 24 25 26 27 


Radio Power 

Figure 8 . Stellar mass versus radio power for RL AGN. The color 
map shows the redshift of the sources. 

fying SFGs as R.L AGN since these objects are clear outliers 
from the radio-FIR correlation, suggesting that there is in¬ 
deed a significant radio emission associated with jets from 
the central black hole. 


8.2 SF in RQ AGNs and SFGs 

Once the contribution from RL AGN is removed, a deep 
radio survey, as the one presented in this work, can be used 
to study the SFG population, its distribution in sSFR, and 
its evolution with redshift. Moreover, we have shown that we 
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can use the radio as an SFR tracer also for the hosts of most 
RQ AGNs, including sources in which the AGN emission 
dominates over the stellar emission in the UV, optical or X- 
ray bands. This allows for a more complete census of SFGs 
and permits to study the AGN contribution in the SFGs 
population. 

Our VLA survey is one of the deepest radio surveys 
available but still it is strongly affected by our flux density 
limit that, as described in Sec[T] allows us to detect the bulk 
of the SFGs population only in the local Universe. There¬ 
fore, it is important to carefully take into account selection 
effects when drawing any conclusion. In this context, a com¬ 
parison with a mass selected catalog ue can be useful. We 
used a mass selected mock catalogue (iBernhard et al.ll2014l f 
of SFGs based o n the empirica l mode l for the SFG popula¬ 
tion d escribed in lSargent et al.l d2012T) and IBethermin et al.l 
(2012). We will now briefly describe the model and the mock 
catalogue before continuing with the analysis of the SF prop¬ 
erties of our radio selected SFGs and RQ AGNs. 


8.2.1 Mock catalogue description 


The two basic ingredients to build the mock catalogue are: 
(i) a description of the observed SFG mass function (MF) 
and its redshift evolution; (ii) a prescription to associate a 
SFR to each source. For the MF we a dopt a para metrization 
similar to th e one used in Sargentet^L_|2012 ) , base d on 
the results of lPeng et al.1 l|201flh and lllbert et, al (j 201 d 'l (see 
Appendix[A]for more details). This allows us to build a mass 
selected SFG catalogue down to a stellar mass of 10' Mg 
(which is smaller than the minimum mass measured in our 
sample) for a mock observation on an area equal to the size 
of our field (0.285 deg 2 ). 

A SFR is then assigned to each source in order to re- 
prod uce the observed MS of SFGs with its redshift evolution 
(e.g. iDaddi et aD 120071 : lElbaz et al.l 120071 : iRodighiero et all 
l201ll : IPeng et al.l 1 201(1) . and the distribution of sSFR at a 
fixed stellar mass. The la tter is modelled as the sum of two 
Gaussian components llSargent et al J[2012h one associated to 
MS galaxies and the second to starbursts. The full width half 
maximum (FWHM) of the two Gaussian, the relative con¬ 
tribution of the two components and i ts evolutio n wit h red - 
shift are free parame t ers. W e refer to ISargent et al] (120121 ') 
and IBethermin et all d2012Tl for a better description of this 
model and to Appendix 0 for an illustration of the model 
parameters adopted in this work. 

In order to directly compare this mock catalogue with 
our radio sample we need to "observe" it with the same 
flux density limit of our VLA survey. To do that, we con¬ 
verted the SFR of the sources in the mock catalogue to the 
corresponding radio power according to equation 0 adding 
a 0.2 dex random scatter, in agreement with what we mea¬ 
sure (see Sec. 0 . We then computed the corresponding radio 
flux density given the source redshift and selected only ob¬ 
jects with radio flux densities above 37 /rJy, the average flux 
density limit of our VLA survey. Finally we added random 
uncertainties with 0.2 dex dispersion on the mass and SFR 
of the mock objects to take into account the uncertainties 
affecting real measurements. 

Sources in the mock catalogue are flagged as starburst 
or MS object depending on which of the two Gaussian com¬ 
ponents they belong to. Recovering this information in real 


observations is not possible for individual sources, but only 
statistically. Therefore, to compare the starburst fraction 
and its dependence on physical parameters (e.g. stellar mass 
or radio power) in the mock catalogue and in our sample, 
we will adopt the same definition of starburst, MS and pas¬ 
sive based on the distance of the source from the MS as 
defined in section [8] Note that even if the mock catalogue 
include in principle only active galaxies, due to the errors we 
added to the mock catalogues quantities and to the wings 
of the MS gaussian distribution, some sources in the mock 
catalogue are "passive" according to our definition. We also 
note that, as the RL AGNs are mainly hosted in passive 
galaxies, the mock catalogue can be used for a comparison 
with the properties of SFGs and RQ AGNs, where the radio 
emission traces the SF activity. 


8.2.2 Starburst contribution 

In our observed VLA sample the large majority of SFGs and 
RQ AGNs are MS objects (~ 74%). The remaining sources 
are mainly starbursts (~ 22%) with a small contribution 
from passive galaxies (~ 4%) (see also Fig. 0. The total 
observed fraction of SB (/sB,obs) is slightly larger in RQ 
AGNs compared to SFGs (25% vs. 20%). We emphasize that 
these fractions refer to our radio-flux limited sample, and 
should not be confused with fractions for a mass selected 
sample. T his result echoes what found for a FIR selected 
sample bv lGruppioni et al.1 l2013l b 

In more detail, we can follow /sB.obs as a function of 
stellar mass for the two classes of sources (middle panels 
in Fig(7J). For both SFGs and RQ AGNs, /sB.obs decreases 
quickly with increasing stellar mass as at low masses we are 
highly biased towards high sSFR due to our flux density 
limit. The behaviour in RQ AGNs is even more extreme 
and the SB galaxies become the majority at masses smaller 
than ~ 1O 1O ' 5 M0. This is mainly due to the fact that we 
are able to detect low mass objects only at low redshift (see 
first panel of Fig. 0, but, since the volume that we probe 
at z< 0.5 is small, the chance of detecting AGNs is very 
low. We are able to sample the MS at low masses only for 
more numerous objects as the SFGs. The effect of the flux 
density limit becomes clearer looking at the right panel of 
Fig. 0 we show the starburst and main sequence galaxies 
relative fraction for the mass selected mock catalogue (light 
colors and empty symbols) and for the same sample after 
applying the flux density cut (dark colors and full symbols). 
In the mass selected case, the SB fraction is constant with 
stellar mass while the "observed" mock catalogue shows the 
same trend as in the real VLA samples. 

More interestingly, we can also follow the evolution of 
/sB.oba, that is equivalent to the relative SB contribution 
to the SFR density, as a function of redshift. Some hydro- 
dynamic simulations suggest an increase of the SB fraction 
out to z~2 due to the increase of the merger-indu ced star 
formation at high redshift (e.g.. lHopkins et al .11 201(1 ). In our 
VLA sample, we observe a fast growth from the local Uni¬ 
verse up to Z'-M followed by a flattening both for RQ AGNs 
and SFGs (see middle panels of Fig. 0. This is not only a 
consequence of our flux density limits but also reflects the 
intrinsic growth of the starburst contribution from the local 
Universe to the peak of the star formation activity. 

Again, the comparison with the mock catalogue allows 
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us to better quantify our results. The evolution of the frac¬ 
tion of starburst in a mass selected sample (/sB.mass), is 
one of the ingredients of the empirical model used to build 
the mock catalogue (see Appendix © and we can therefore 
use our observation to constrain the model parameters. We 
fixed the loca l value of /sB.mass to 1.2% as computed by 
ISargent et al . 1 (f20T3> as our small field does not allow us to 
have enough volume to properly constrain it, while varying 
the slope of the growth and the redshift above which the 
evolution stops. We find that a growth as (1 + z ) 2 up to z=l 
of /sB.maas and constant thereafter is required to match the 
total /sB.oba in our combined sample of SFGs and RQ AGNs 
and to reproduce its behaviour as a function of redshift (see 
Fig. [6]). 

In summary, the "observed" mock sample reproduces 
well the behaviours of the relative fraction of SB, MS and 
passive galaxies in our SFGs, both as a function of redshift 
and stellar mass (see the last two panels of Fig. [5]and Fig. [7] 
respectively). We note that also when combining SFGs and 
RQ AGNs, that are the minority in number, the fraction 
of SB as a function of stellar mass and redshift is consis¬ 
tent with the model within the uncertainties. Also the to¬ 
tal /sB.obs is the same in the VLA and mock flux limited 
samples (~ 20%). Therefore, we conclude that the empirical 
model of lSargent et al. lfloll) is suitable to describe the star 
forming galaxy populations with only a little tuning of the 
model parameters. 

Therefore, we could use the model prediction to study 
the behaviour of a mass selected sample rather than a flux 
limited ones. In the last panels of both Fig. [S]and[T]we show 
in light colors and empty symbols the fraction of SB and MS 
galaxies for the mass selected catalogue. As expected, the 
"real" /sB.mass is much lower compared to the one observed 
in a flux limited sample and SB sources contribute to less 
than 10% to the total SFR density. This resu lt is al so in 
agreement wi th the results o f iRod ighiero et al.l (1201 il l and 
the models of lHopkins et al.l (l2010h . 


8.3 AGN content as a function of sSFR 

In Fig. [9] we plot the relative fraction of the dif¬ 
ferent classes of sources as a function of the dis¬ 
tance with respect to the MS (i.e., A log(sSFR)MS = 
log[sSFR(galaxy)/sSFR,Ms(M s ta.r,z)]). Below the MS 
(left side of Fig. 0 the population is dominated by RL 
AGNs. This is due to two main reasons. The first is re¬ 
lated to selection effects since we are sensitive to low sSFR 
only for sources whose radio emission is enhanced by AGN 
activity, i.e., RL AGNs. On the other hand, the probabil¬ 
ity of hosting a RL AGN has b een observed to be higher in 
massive passive galaxies (e.g., iBest et akl 120051 ) even if its 
physical reason is still unclear. 

On the right side of the plot, i.e. above the MS, we ob¬ 
serve a rising of the RQ AGNs fraction from ~25 to 40%. To 
investigate if this rising can be due to some AGN selection 
effects we considered two SFG sub-samples matched in red¬ 
shift and mass with the RQ AGN distributions. The trend 
remains in both cases. We also checked if it can be explained 
by some AGN contribution to the FIR luminosity in two 
ways: (i) using the SFR derived from the radio power rather 
than SFRfib, and (ii) using the SFRfui corrected for the 
AGN contribution (as describe in Sec l6.ll) to compute the 



Figure 9. Relative fraction of sources classes as a function of 
A log(sSFR) MS . 


Alog(sSFR)MS- Again the fraction of RQ AGNs increases 
as the Alog(sSFR)MS in creases. This resu l t is also in agree¬ 
ment with what found bv lGruppioni et alJ ( ; 2(~1 .' I). who con¬ 
sidered an FIR. selected sample and found that starburst 
sources are dominated by galaxies with AGN-type SEDs. 
We will further expand on this result in a coming paper 
(P. Padovani et al. in prep). The most extreme SF activity 
seems therefore often associated with an active phase of the 
black hole. Such scenar io has been predicted by simulation of 
merging systems (e.g., iDi Matteo et al.ll2005l : Hopkins et alJ 
120081) : the instabilities generated by such a violent event trig¬ 
ger both a burst in the SF and feed efficiently the central 
BH. 


As discussed in the previous section, sources with SB 
activity represent a small fraction of the SFGs population. 
Similarly, the large majority of our RQ AGNs (75%) are 
hosted in MS or passive galaxies. This fraction only repre¬ 
sents a lower limit since our observation do not fully sample 
the MS especially at high redshift. Hence, the most common 
triggering mechanisms of the black hole activity is proba¬ 
bly related to steady process like gas inflows from the in¬ 
ter galactic medium rather than extreme events like m ajor 
mergers (e.g., iMainieri et al.ll201ll : iRosario et al.ll2012l h 


9 DISCUSSION 

We find that radio power is a tracer for the SFR as good as 
the FIR luminosity not only for SFGs, but also RQ AGNs. 
Hence, our finding suggests that radio emission in RQ AGNs 
is mainly due to SF. As already mentioned in Section[ 6 ] even 
if some level of correlation between the radio and FIR. lu¬ 
minosity was expected due to the fact that RQ AGNs lie 
in the same locus as SFGs in the 524068 -redshift plane, the 
small scatter we find, of only 0.23 dex, was not foreseeable. 
Indeed, in the MIR, that we use to separate the RL AGNs 
(see section 12721) . the AGN contribution can be strong lead¬ 
ing to a scatter in 52406 s values for RQ AGNs two times 
larger compared to SFGs. Nevertheless, we have shown that 
using only two bands, the radio (1.4 GHz) and the MIR 
band (24 pm), we can efficiently and with few outliers sep¬ 
arate sources powered by the two different radio emission 
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mechanisms, namely jets and SF. Since the UV-optical and 
also the MIR emission can be heavily contaminated by the 
AGN emission in powerful AGN, like most of our RQ ob¬ 
jects, the radio power could provide a better estimate of the 
SF in their host galaxy compared to e.g., UV-based tracers. 
This result can be particularly useful especially for deep ra¬ 
dio surveys without enough FIR. coverage. In addition, being 
radio frequencies almost unaffected by dust extinction, it is 
suitable for both type I and type II AGNs. 

Hypothetically, the good agreement between SFR T and 
SFRfc for RQ AGNs could be the result of a conspiracy 
where the contribution to the radio luminosity from tiny 
jets is exactly compensated by an extra contribution in the 
FIR. This will boost both the SFR estimates by the same 
amounts. With the current resolution of the VLA data we 
cannot exclu de the presence of jets in the center of these 
AGNs fe.g.. [Giroletti fe Panessa 120091 1 and higher resolu¬ 
tion radio observations with the VLBI would be needed to 
spatially resolve them. 

9.1 Comparison with a FIR selected sample 

We want to investigate the possibility of using deep radio 
surveys as a useful tool to study the cosmic SF history, al¬ 
ternative or complementary to FIR surveys. We have there¬ 
fore compared our sample with the Herschel /PACS detected 
sample. 

As described in section 12731 only ~60% of our VLA sam¬ 
ple has a counterpart in the Herschel catalogue of the E- 
CDFS. Excluding the RL AGNs, for which we expected a 
low detection rate as most of them are hosted in passive 
galaxies (but see Sec. 18.111 . the fraction increases to ~ 70%. 
As for the 30% of unmatched radio sources, this is at least 
partly due to the dispersion of the RFC: a radio and a FIR 
flux density limited sample would tend to be biased towards 
sources on opposite sides of the relation. To test this hypoth¬ 
esis, we made use of the mass selected mock catalogue de¬ 
scribed in Sec. 18. 2. II We "observed" it with a FIR-luminosity 
limit equivalent to our 37/rJy flux density limit computing 
the minimum SFR corresponding to the radio flux density 
limit at each source redshift and selecting only the objects 
with SFR above this threshold. The resulting FIR-selected 
catalogue contains ~780 sources, roughly the same number 
of objects as in the radio selected mock catalogue described 
in Section 18.2.11 In Fig. [TO] we show the RFC for the FIR 
selected mock catalogue in cyan, the radio selected in ma¬ 
genta, and for the matched sources, in blue. The latter are 
only 60% of the total number of objects in each catalogue, a 
fraction similar to the one obtained matching the real VLA 
and Herschel samples. The offset from the 1:1 correlation is 
~ 0.15 dex above and below for the FIR and radio selected 
samples, respectively. This plot explains also why most of 
the VLA sources without PACS detection are just below 
the correlation found for the VLA sources with Herschel 
counterpart in Fig. [3] being radio-detected only they have a 
slight radio excess. 

On the other hand, about 40% of the PACS detected 
sources have a counterpart in the VLA image. This means 
that the FIR observations reach a lower flux density limit, 
also thanks to the much deeper maps in the central GOODS 
held compared to the outskirts ( see Sec. I m We have c om- 
pared the SFRfir computed in lGruppioni et all (120131') for 



SFR radio (M ( ~, yr ’) 


Figure 10. SFR. derived from the FIR luminosity versus the SFR 
from the radio luminosity for the sources in the FIR selected mock 
catalogue (cyan), the radio selected mock catalogue (magenta) 
and the matched sample (blue). The dashed lines are the linear 
fit to the data. 

the sources in common and we find very good agreement 
with our estimates. Also the stellar masses are consistent; 
we find a larger scatter but no systematic shifts. This re¬ 
flects the fact that the SFR is well constrained when FIR 
photometry is available, while on the stellar masses there 
are larger uncertainties, especially in AGN host galaxies, de¬ 
pending on the fittin g techniq ue and especially on the model 
library adopted (e.g. filbert et al.ll20ld ). 

In Fig. 1111 we show the VLA (RQ AGNs and SFGs 
only) and the Herschel samples on the SFR-stellar mass 
plane with the sources coloured according to their redshift. 
The two sets of data show the same trends and the same 
evolution with redshift. We also find a similar fraction of 
SB (20%) in the PACS sample compared to the radio one. 
Being a bit deeper, especially in the GOODS field, the PACS 
observations probe better the small stellar mass end and are 
able to detect more low SFR sources at low redshift. Already 
ongoing deeper radio surveys, as well as those planned with 
future faciliti es, w ill allow to p robe these small masses also 
in the radio IlNorris et al.ll2013l) . Indeed, going down to the 
nanojy sensitivity level, they will be able to probe the bulk 
of the SFGs population. An illustration of that is given in 
Fig. [T2] where we plot the fractions of MS and SB galaxies 
as a function of the radio flux density in /j ,Jy for the mass 
selected mock catalogue. Below about 10 /rJy the number 
density of SB remains nearly constant around 8%, while at 
the flux density limit of our survey (vertical line) we are 
biased towards higher sSFRs. 

9.2 What changes using a different SFR tracer? 

As discussed in the previous section, due to the RFC disper¬ 
sion, a sample selected with a radio flux density limit would 
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Figure 11. SFR versus stellar mass for the VLA sample (left) and for the PACS sample (right). 


Mass selected mock catalog 



Figure 12. Fraction of MS and starburst sources in the mass 
selected mock catalogue as a function of the radio flux density. 


be slightly biased towards those objects on the right side 
of the relation, i.e. with a radio power larger than the one 
given by the relation for a given FIR luminosity. 

That means that, especially for the sources without 
Herschel counterpart, the SFR. derived from the two dif¬ 
ferent tracers can be slightly different. As quantities like the 
fraction of SB are extremely sensitive to these changes due 
to small number statistics, we investigated how the choice 
of a different SFR tracer, namely the radio power, can af¬ 
fect our conclusions. Of course that cannot be done for RL 
AGNs as in these objects the radio power is not tracing the 
SF activity in the host galaxy but it is highly contaminated 
by the jets emission. For RQ AGNs and SFGs instead we 
repeated the same analysis as in Sec. [8] using the SFR T to 
compute the A log(sSFR)MS and therefore to determine the 
level of their SF activity. 

The total fraction of SB in this radio flux-limited sample 
increases to ~309<@ but the trends as a function for example 
of stellar masses remain the same (see Fig. ED). Also in the 
radio flux density limited mock catalogue, if we compute the 
A log(sSFR)MS from their radio power we obtain a /sB.obs 
of 30%, in agreement with our data. 


4 Note that, in contrast, the fraction of passive goes almost to 
zero. 


The larger differences appear at z > 1, where the frac¬ 
tion of starbursts is up to a factor of two larger using the 
SFRi, but this effect is mainly due to the sources without 
PACS detection. Indeed, as shown in Fig. [4] the SFR ex¬ 
trapolated from the optical-to-MIR photometry tends to be 
underestimated at high redshift. Hence, the differences be¬ 
tween the two SFR tracers appear larger. As a further con¬ 
sequence, the /sB,obs at high redshift shown in the central 
panels of Fig.[|J]are slightly underestimated as also suggested 
by the predictions of the empirical model. 


9.3 Lx — Lfir relation 


The hard band X-ray luminosi ty (Lx) is considere d another 
SFR tracer up to high redshift llNorman et al.ll2004l . In star- 
burst galaxies the 2-10 keV rest-frame emission is dominated 
by high-mass X-ray binaries (HMXB). Being the compan¬ 
ion of the accreting object a short-lived high mass star, the 
Lx is linked to the recent star formation activity. At low 
star formation rates, however, there is contamination from 
low-mass X-ray binaries (LMXB) which do not trace the im¬ 
mediate SFR due to their long evolution lifetimes. The con¬ 
version factor between Lx and the SFR has been calibrated 
for local samples (e .g. jRanalli et al.ll2003l : IPersic et alj|20()3 : 
iLehmer et alj[2012h . It has been argued that the hard band 
X-ray luminosity can be used as an unbiased SFR indica¬ 
tor up to high redshifts, despite the contamination from 
obscured AGN which may be unnoticed in distant, faint 
sources._ 

Recently. I Vattakunnel et al.l (l2012f) used the deep X-ray 
Chandra data (total integ ration time of 4 Ms, see lXue et al.l 
|201lt_ LehmeretalJ 12012 ) and the deep VLA observation 
I Miller et al ■2S" to explore the Lx-LV relation in galaxies 
dominated purely by star formation processes in both bands. 
Due to the strict selection process, aimed at discarding any 
significant AGN contamination, they obtained a sample of 
43 SFGs up to z ~ 1.2. They found a clear linear correlation 
between radio and X-ray luminosity in SFG over three or¬ 
ders of magnitude in this redshift range, consistent with that 
measured locally. They also measured a significant scatter 
of the order of ~ 0.4 dex, higher than that observed at low 
redshift, implying an intrinsic scatter component, or some 
residual AGN contamination. 

The information on the IR properties of the E-CDFS 
sources presented in this work, allow us to further explore 
the relation between Lx and SFR in this sub-sample of 
sources. As a first step, we consider only 35 sources, out of 
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RQ AGN SFG Mock catalog 





Figure 13. Fraction of starbust (magenta circles), MS galaxies (green diamonds), and passive galaxies (red squares) derived using the 
SFR r , versus stellar mass. The panels refer from left to right to RQ AGNs, SFGs and mock catalogue sources. 


43, which have a > 5cr radio detection, and plot the Afir- 
P r relation. We find that the bulk of the sources follow the 
radio-infrared correlation expected for star forming galaxies 
(see Fig. 1141 left panel), confirming that the selection was 
effective in selecting sources powered by SF processes. How¬ 
ever, we find 4 clear outliers at low IR luminosities, that 
indeed we re classified as RL AGNs according to their q^Aoba 
value (see lBonzini et ahll2013l . for details). For the remaining 
31 sources, we plot the Ax-Afir relation (see Fig. 1141 right 
panel). 

We also measure the observed scatter in the Ax-Afir 
and Afir-AV relations finding that in the first case the dis¬ 
persion is ~0.29 dex, about two times larger than the one 
between Afir-Rc (0.17 dex). Even after removing the statis¬ 
tical error associated to the X-ray and radio luminosity, the 
intrinsic scatter of the Ax-Afir relation remains two times 
larger than the other. Given the robust selection, we can 
safely assume negligible AGN contamination, and therefore 
interpret this scatter as an intrinsic scatter associated to the 
SFR tracer. Therefore, for the first time we show that the 
X-ray luminosity, as a star formation tracer, is noisier than 
radio and IR luminosities. 

The scatter can be due to an increasing X-ray emis¬ 
sion component proportional to the stellar mass and not to 
the instantaneous SFR, analogous to what has been found 
in a s ample of luminous infrared galaxies by iLehmer et al.l 
d2Qloh . Alternatively, the scatter may be due to the hard 
component of a strongly absorbed AGN component. 


10 SUMMARY 

In this work, we investigated the star formation properties 
of the faint radio population as detected by one of the deep¬ 
est 1.4 GHz survey up-to-date conducted with the VLA 
in the Extended Chandra Deep Field South . This stu dy 
buil ds upon the results presented in iBonzini et al.l (120121 1 
and lBonzini et akl ll2013l ) where we have exploited the wealth 
of multi-wavelength data available in this field to identify 
the AGNs, further divide them into radio-loud and radio¬ 
quiet, and to characterize the properties of the radio selected 
galaxies (e.g. redshift, stellar mass). The main results of this 
paper are the following: 

• We have derived the FIR luminosity, fitting the UV-to- 
FIR SED of our VLA sources; radio selected SFGs follow 


the RFC up to z~3 with a nearly constant dispersion of 0.2 
dex (Sec.[4]) 

• Comparing the SFR derived from the FIR luminosity 
and the radio power (Tabic [Bill , we show that the two SFR 
tracers are equivalently good not only in non-active SFGs 
but also for the host galaxies of RQ AGNs. This implies that 
the main contribution to the radio emission in RQ AGNs (at 
least at 2 ~ 1.5 — 2) is associated with SF activity in the 
host rather than with radio jets (if present) powered by the 
black hole (Sec. 0 

• As the SFR in SFGs correlates with their stellar mass, 
we made use of the sSFR to determine the SF activity 
level in our sample (Sec. [8j. The data are consistent with 
models that predict that the SB galaxies represent a small 
fraction (^8%) of a mass selected SFGs p opulation (e.g., 
ISargent et all 120121 : iBethermin et all l2012h . We have also 
discussed the impact of our flux density limit on a mass 
selected mock catalogue (Sec. 18.21) . 

• The majority of our RL AGNs are hosted in passive 
galaxy, but we detect significant SF activity in ~40% of the 
RL AGNs hosts. This suggests that, at least for low radio 
power, the presence of radio jets does not always prevent 
active star-formation (Sec. 18. ill . 

• We find hints of a higher fraction of AGNs in the most 
extreme SB galaxies (Sec. [H751) . This is in agreement with 
scenarios where both the intense SF activity and the efficient 
accretion on t he black hole are triggered by gas rich m ajor 
mergers (e.g., IChen et al.ll2013l : ISvmeonidis et al.ll2013l ). 

• The vast majority (~75%) of our RQ AGNs lie along 
the main sequence suggesting that the bulk of the black hole 
activity is associated with secular processes. 


Finally, we have shown that deep radio continuum surveys 
are a powerful tool to investigate the star formation history 
up to high redshift. The current and up-coming radio facil¬ 
ities like the JVLA and the Square Kilometre Array (SKA) 
pathfinders will be able to ob serve with nanojy (lcr) sen¬ 
sitivity large areas of the sky ([Norris et al.ll2013| ) hence de¬ 
tecting the bulk of the SF population. Therefore, our finding 
that the main contribution of radio emission in RQ AGNs is 
due to star formation in their host galaxy opens the possi¬ 
bility to use radio emission to estimate the SFR even in the 
host galaxy of bright quasars. 
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Figure 14. Left: P T — Ljr relat ion for the sub -sample of SFGs with X-ray counterpart selected in IVattakunnel et al.1 (|2012h . The four 
sources classified as RL AGNs in Bonzin i et all d2013T) are plotted as red squares. Right: Lx — Ljr relation for the same subsample after 
removing the outliers in the P r — L \ |< relation. 
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APPENDIX A: CHOICE OF THE MODEL 
PARAMETER 

We describe here in more de tails the empirical m odel used 
to build the mock catalogue ijBernhard et al.ll2014i l and dis¬ 
cuss the set of model parameters adopted. The physical mo¬ 
tivation fo r the model and the und erlying equations are de - 
scribed in lSargent et al.l (120121 1 an d lBeithermin et ah] ll2012lh 

The motivation for thi s Appendix is that the original 
set of parameters adopted in I Bethermin et needed 

to be modified in order to reproduce our observations. In¬ 
deed, the corresponding mock catalogue, "observed" with 
the same flux density limit of our survey (see Sec. 18.2.11) . 
while reproducing approximatively the radio counts, does 
not well reproduce some host galaxy properties of the sam¬ 
ple; the mock catalogue has a mass distribution that peaks 
at higher masses compared to the distribution for our VLA 
sources. As a consequence, also the SFR properties are not 
well reproduced in the mock catalogue having on average 
lower SFR and therefore a significantly lower fraction of star- 
burst galaxies compared to the data. 

Therefore, we investigated the possible reasons for this 
discrepancy and looked for a set of parameters that is able 
to reproduce our observations in terms of mass and sSFR 
distribution together with the radio counts. 

Al SFGs mass function 

The mass distribution of the mock catalogue objects is set 
by the SFG mass function (MF)._The MF adopted in the 
model is based on the fits b vIPeng et alj (l2010h of the SFGs 
MF presented in (Illbert et al .11 2Old ). It is described as a sin¬ 
gle Schechter function with characteristic mass (Mb) and 
faint-end slope that are redshift invariant and by a constant 
characteristic density up to z~ l followed by a decline as 
(1 — z) 0 ' 46 Isee lBethermin et al.l (l2012l l for details). Our ra¬ 
dio observations are relatively shallow compared to the K- 
band data used to compute the MF in Illbert et alj (l2010h . 
consequently we are more sensitive to the high mass end of 
the distribution. Therefore, we are particularly sensitive to 
the value of Mb, that set the position of the MF break fol¬ 
lowed by the exponential cut off. Since our sample is small 
and highly incomplete at low masses, we cannot perform 
a real fit of the MF. We chose to keep all the other pa¬ 
rameter describing the MF fixed. As already mentioned our 
observations have a mass distribution that peaks at lower 
masses compared to the corresponding mock catalogue. To 
reconcile the model with our observ ations we assume a Mb o f 
lO n M 0 rather than 10 11,2 Mq as in lBethermin et all (|2012|J . 
Adopting this value the mass distributions of the mock cat¬ 
alogue and of our data are consistent as confirmed by a 
Kolmorogov-Smirnov (KS)-test (Prob> 0.99). The possible 
reasons for this shift are twofold: on one hand, stellar mass 
measurements have large uncertainties and can be up to 0.3 
dex systematically different depending on th e ste llar popula¬ 
tion synthesis model adopted (see e.g. Ilbert et al.ll201(l ). To 
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check this hypothesis we computed the stellar mass us in g the 
same method adopted for the VLA source s j Bonzini et ah 
1201 .'l l for a sub-sample of the lllbert et, al.1 (|201f)h sourced 
and com pared the mass m easurements. We indeed noted 
that the lllbert et al.l (|20ld ) stellar masses are on average 
larger than the one obtained with our method. On the other 
hand, there is no one-to-one correspondence between the 
color-c olor based m e thod used to identify the SFG popula¬ 
tion in lllbert et al.l (l20ld l and our scheme; it is therefore 
possib le that some sources excluded from the lllbert et all 
(12010) SFG sample because of their redder colours have been 
considered SFGs according to our classification. Finally, we 
note that the simple prescription adopted in the empirical 
model and the set of parameters chosen is only an approxi¬ 
mation for the S FGs MF while a more complex description 
as e.g. in lllbert et al.l (120131 ) would be needed. 


A2 Main sequence and its redshift evolution 

T he two mod es of SF are described by two Gaussians 
llSargent et al .1 12()1 2l i. We assume a FWHM of 0.2 dex for 
the MS galaxies as measured in lRodighiero et al.l (1201 il l. For 
the second Gaussian w e use the same parameter adopted in 
ijBethermin et al1l2012l ). i.e. a FHWM of 0.2 dex, a displace¬ 
ment from the MS peak of 0.6 dex. The relative fraction of 
SB evolves with redshift as: 


fsB = 0.012 x (1 + z) 2 for z<l (Al) 

The local fraction is taken from 1 Bet hermin et al.1 ( 2012 1 but 
we set a steeper growth of the starburst fraction up to red¬ 
shift 1 as suggested by our observations (see Sec. IQt) . 

The normalization of the MS is another critical parame¬ 
ter in determining both the number counts and the observed 
SB fraction. Keeping fixed the slope of the MS to the value 
measured in lRodig hiero et alJ (1201 ill and its evolution with 
redshift to ( 1+z ) 2 8 (jSargent et al.ll2012l ). we looked for the 
best value for the MS normalization. We adopt a normal¬ 
ization of the MS ( loq sSFR(z = 0, M = 10 11 Mq) = 10.08) 
that is lower than in I Bet hermin et al . (l201lj) b ut it is in 
better agreement with what was found in Rodighiero et all 
(12011 1 and our own measurement. It is also c onsisten t with 
previous results in the literature up to 2 : ~2 dNoeske et al.1 
l2()0?l : [E~1baz et aI]|2007l : rDaddi et al .112007?) . Recent work sug¬ 
gests that above 2 ~2 the evolution of t he MS is sha l¬ 
lower than (1 + z) 2 ' 8 as assumed by th e Bemharde^lJ 
(|2014|) model (e.g. Weinmann et al.ll201ll : IStark et alJ 20131 : 
iGonzalez et al.l [201 dl l. As a consequence, the MS shown in 
the last panel of Fig. [5] is probably too high. However, 
only about 5% of our sources fall in this high redshift bin 
(2.5< * <4) and therefore the simplification adopted in the 
model does not affect significantly our conclusions. We stress 
that the main goal of the comparison of our data with the 
mock catalogue was to better quantify the selection effects in 
flux density limited samples, to check the qualitat ive agree¬ 
ment with the two Gaussian SF mode model dSargent et al.1 
121)1 J i and not to obtain a fit of the MS as a function of 
redshift. 


A3 Comparison with VLA observations 

With the few changes in the model parameters described 
above, we are able to reproduce both the radio counts and 
the physical properties of our radio sources with the mock 
catalogue. Indeed the mass selected mock catalogue, "ob¬ 
served" with the same flux density limit of our VLA survey 
(see Sec l8.2.lll , on a mock field with the same area of our ob¬ 
servation, contains about 790 sources. The exact number of 
sources varies of some units for different runs of the model 
since we add random uncertainties on the mock galaxies 
properties and a random dispersion for the RFC. The VLA 
sample considered in this work contains 779 sources. Consid¬ 
ering only the RQ AGNs and the SFGs, where the main con¬ 
tribution to the radio flux is due to SF, but c orrecting for the 
not p erfect uniformity of the sky coverage llPadovani et al.1 
I2015T) . the number of RQ AGNs and SFGs detectable at the 
37/rJy flux density limit in our field is 784. Considering that 
some of the RL AGNs, according to our analysis, contribute 
to the SFG population, we conclude that the number counts 
predicted by the empirical model are in good agreement with 
the observed one. A complete analysis of the impact of the 
various model parameters on the mock catalogue physical 
characteristics is beyond the scope of this paper. The goal of 
this investigation was to check the consistency between our 
ob servation and the model for the SF population proposed 
in ISargent et al.l d2012l ) and to better control the effects of 
the flux density limit on our results. 


APPENDIX B: STELLAR MASSES AND SFRS 
CATALOGUE OF VLA SOURCES 

We make publicly available the physical properties derived 
for our radio sample and used in this work. They are sum¬ 
marized in table [nu The catalogue columns are organized 
as follows: 

• (1) Identification number of the radio source (RID). 

• (2) Source classification. 

• (3) Source activity: "SB" for starburst galaxies, "MS" 
for main sequence galaxies, and "P" for passive galaxies ac¬ 
cording to the definition given in section [8] 

• (4) Source redshift. 

• (5) Stellar mass. 

• (6) SFR derived from the radio power. 

• (7) SFR derived from the FIR luminosity. 

• (8) distance with respect to the main sequence in the 
SFR — M a tar plane (Sec. 18.31) 

• (9) PACS flag; 0= only upper limits in the PACS bands, 
1= at least one PACS detection 


5 A sample of X-ray selected AGNs for which we have optical-to- 
24/rm photometry and redshift. 
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Table Bl. Star formation properties of the VLA sources.(The full Table is available in the on-line journal. A portion in shown here for 
guidance regarding its form and content.) 


(1) 

id 

(2) 

class 

(3) 

type 

(4) 

z 

(5) 

Mstar 

[log{M Q )\ 

(6) 

SFRr 

\log(MQyr~ 1 )] 

(7) 

SFR Fm 

[log^MQyr- 1 )} 

(8) 

sSFRex 

(9) 

PACS? 

711 

RL AGN 

P 

1.89 

10.774=0.18 

277.94=77.6 

44=0.1 

-1.30 

0 

712 

SFG 

MS 

0.56 

10.794=0.05 

15.44=4.7 

214=1.1 

0.07 

1 

713 

RQ AGN 

MS 

0.49 

10.604=0.06 

17.24=3.3 

104=0.8 

-0.01 

1 

714 

SFG 

MS 

0.77 

10.304=0.02 

41.04=9.3 

304=4.0 

0.47 

1 

715 

SFG 

SB 

0.18 

9.954=0.07 

13.34=0.8 

114=0.2 

0.81 

1 

716 

RQ AGN 

MS 

1.16 

10.784=0.10 

80.24=25.1 

794=20.9 

0.27 

1 

717 

SFG 

MS 

0.52 

10.134=0.02 

17.74=3.8 

124=0.9 

0.41 

1 

718 

SFG 

SB 

2.46 

10.524=0.09 

478.74=141.2 

7034=86.2 

0.84 

1 

719 

RL AGN 

MS 

1.03 

10.844=0.19 

38394.14=20.2 

194=0.2 

-0.31 

0 

720 

SFG 

MS 

0.25 

10.594=0.02 

16.54=2.4 

104=0.5 

0.19 

1 
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